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Abstract 

We study solutions of the relativistic hydrodynamical equations, 
I which describe spherical or cylindrical expansion of ideal fluid. We 

derived approximate solutions involving two arbitrary functions, which 
describe asymptotic behavior of expanding fireballs in ultra-relativistic 
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O limit. In case of a linear equation of state p(e) = ke — c\, (0 < k < 1) 

| we show that the solution may be represented in form of an asymptotic 

series in negative powers of radial variable; recurrence relations for 
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the coefficients are obtained. This representation is effective, if k > 
1/(2J + 1) (J = 2 for spherical expansion and J = 1 for the cylindrical 
one); in this case the approximate solutions have a wave-like behavior. 



^ ■ 1 Introduction 

The problem of the relativistic ideal fluid expansion arise, e.g., in the theory 
of gamma-ray bursts, that represent the most powerful explosions in the 
Universe p. The other applications are due to hydrodynamical theory of 
multiparticle production j2] . Here the relativistic ideal fluid model is utilized 
to describe the behavior of relativistic fireballs. Detailed investigation of 
this problem involves complicated 3-D numerical simulations. Nevertheless, 
considerable information may be obtained from analytical studies, which deal 
typically with one- dimensional flows having spherical, cylindrical and plane 
symmetry jHj, [2]. One of such well-known results is due to Blandford & 
McKee j3j, who studied self-similar ultra-relativistic gas expansion following 
a strong shock. Blandford & McKee solution is widely used for interpretation 



1 



of the gamma-ray bursts and their afterglows pQ. Certain drawbacks of this 
solution have been pointed out in [5]. 

It should be noted, that investigation of ultra-relativistic asymptotics of 
ideal flows is rather specific; it requires some care to provide approximation 
correctness. For example, if we leave only the main terms in the hydro dy- 
namical equations in ultra-relativistic limit, these become degenerate (see 
Sect. 2). To be sure about approximation involved, one must study the 
higher order corrections in the equations. Evidently, it is desirable to have 
a sufficiently general form of the solutions that contain necessary number of 
arbitrary functions. This is important, e.g., in order to describe models of 
expanding fluid with different radial density profiles. 

In this connection we propose an approximation method to construct 
general solutions that describe ultra-relativistic radial expansion of ideal fluid 
with necessary degree of accuracy. Most detailed result is obtained in case of 
a linear equation of state (EOS); in this case we represent the solution in form 
of asymptotic series (Sect. 3), which is workable at least for a sufficiently 
stiff EOS. These results are compared in Sect. 4 with a self-similar solutions. 



2 Main equations 

The equations of motion of ideal relativistic fluid follow from the conservation 
laws (see, e.g., (2)13) 



fJ,U. 



U 



is the fluid 



d v T» v = 0, 

for the energy-momentum tensor — (p + e)u^u u — pg 
four -velocity, g^ v = diag(l, —1, —1, —1), p is the pressure, e is the invariant 
energy density, the light speed c = 1. We confine ourselves to EOS p = p{e). 
Some kind of the continuity equation must be addedln in a more general case 
of two-parametric EOS. 

In case of cylindrical (J = 1) and spherical (J = 2) symmetry of the flow 
we deal with only two independent variables (time t and radial variable r). 
The equations take on the form [3], [2]: 
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where u r is the radial component of the four-velocity, the other spatial 
four-velocity components are equal to zero. 

In ultra- relativistic limit ti r « u' » 1. But if we leave in Eqs.Q and 
(j2J) only the higher order terms in u r , then both from ((TJ) and (J2J) we have 
the same equation 
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that is we have certain degeneracy. We can determine {u l ) 2 {e + p) from 
this equation, but we cannot determine the velocity and the energy density 
separately. This shows that we cannot neglect the terms, looking at first 
sight as less important, without more detailed consideration. 

3 First approximation 



Further we use the substitution u f = u + l/(4w), u t = u — l/(4tt). Then Eqs. 
fl]} and © yield: 
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Consider the sum and difference of Eqs. Q and (@J) in new variables 
a = t — r, P = t + r and . We have 



d_ _ d fe + p\ J (e + p) f J_ 
dp 1 P> da V Au 2 ) p - a \ Au 2 



d_ 

dp 



u 2 (e + p) 



d (e — p 

+ ~z 



da V 4 



J (e+p) 
P — a 



1 



u 



(5) 
(6) 



These are still exact equations. 
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In the ultra- relativistic approximation u » 1 we leave the terms ~ u 2 
and zero orders in this value. Neglecting the terms containing ~ 1/w 2 , we 
have from © 
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(3 — a 



(7) 



This equation can be solved, if we know the state equation p = p(e). In this 
case Eq. (J7J) yields 



^1 — dp/de] 



de J [e + p (e) 



Let by definition 



$(£:) = exp 



d(3 13 -a 

1 — dp/de 
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de 



e + p{e) 

Then the solution of Eq.(jHJ) may be represented as 

where f(u) is an arbitrary function. The validity of this approximate solution 
depends on smallness of the terms in Eq. (jSJ) that have been neglected, i.e. 
the approximation is valid if 



As we shall see in the next section, this condition is satisfied, e.g., in case 
of a linear EOS p(e) — ke — c x , if 1 > k > 1/(2 J + 1). 
Eq.© can be written as 



((3-a)- J ^[u 2 (f3-a) J (e + p) 



J(e + p) d e — p 
4(/3 - a) ~ da 4 



Using e from Q we obtain u from 

u 2 (p + e)((3-a) J 



g{a) + / dp 
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where g(a) is an arbitrary function. This enables us to find e and u separately. 

If the integral in the right-hand side of (jl(J|) converges for (3 — > oo and 
a remains bounded, then for sufficiently large (3 ~ t we obtain that u 2 (e + 
p)(P — a) J is a function only of a. We shall see below that this is just the 
case of the linear EOS with k > 1/(2 J + 1). In this case we have that the 
energy profile of the ultra-relativistic expanding shell dE/da ~ (e + p)u 2 r J 
remains almost constant for large values of t. 

4 Higher orders of asymptotic series for the 
solutions 

Consider the higher orders of approximation in case of a linear equation of 
state p(e) = ks — ci, where k = Cq < 1, Co is a speed of sound. We are 
looking for the asymptotic representation of the solution, which is workable 
for bounded a = t — r and t — > oo; this means also (3 ~ r —>■ oo. 
We put 

— + Cl /(l + «) ~- J ( 1 + K ) 



(f3-a) a ,y n 1-k 
u 2 (e + p) 
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This substitution reduces the system ©, © to the form 

dip (1 + /Q 2 1 d fip 2 \ 
d(3 + 4 (1-k) (/5- a ) a - J ^ J ~ 

Jk(1 + k) 2 1 <^ 2 
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We shall look for solution of exact equations (jll)) . (jl2j) in the form of 
asymptotic series (m, n — 0, 1, 2, ...) 



*( a >fl = Z j^y!L > (13) 
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(14) 



where 7 = a — J — 1 . 

If we are looking for an asymptotic solution for r — > 00, this representation 
is effective for 7 > 0, that is for k > 1/(2 J + 1). 

We also introduce 



E 



^ ^ (/? ~ a) 



the coefficients of these series may be expressed be means of <p n ,mi ip n ,m ■ 
Moreover, Xn,m depends only upon cp n ' )m ', Tp n ',m> with n' < n, m! < m. 

We now substitute the above representations for ip , tp, x into Eqs.(|TT|). 
(|T2*|) to find recurrence relations for y? n>m , ^n.m- 

For n > 1 , m > 1 we have 
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4 da 
For n = 0, m = 1, 2, ... 



whence also xo,m = 0. 
For m = 0, n = 1, 2, ... 



(771 + m) ^n,m = 

H -7— [7 (n - 1) +m - 1] (p n -i,m-i- (16) 



¥0,m = 4>0.m = 0, (17) 



(1 + K) 2 dXn-1,0 / 10 x 

= 4(1-.) da (18) 

and 

7^n,0 = - (19) 
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Using (|15 p .([16 p . (|T7jl . (jTBJ) we can express all the coefficients of the series 
(HHJ) , (HU) by means of V'oo(a) and <^oo(«) (V'oo 7^ 0). 

Here we write expressions for some lowest order coefficients 

_ (1 + k) 2 dxo,o , _ 1 - k dipo fi 
< ^ 1,0 47 (1 — k) da 1 1,0 47 da ' 

{l + K) 2 dxi,o _ (f h0 (f 0fi Vi,oVo,o 

^2,0 — 7r-q r j — > Xi,o — 1 — ; 72 > 

87(1 - k) da ipQ >0 Vo,o 

(l + K) 2 rf 2 Xo ,o J/t(l + «:) 2 

^2,0 — — l~^2~i ^1,1 — 77 7T7 —TXo.O^ 

2Y da 2 (1 - k) (7 + 1) 

ip ltTn = 0,m = 1,2, 



5 Self -similar solutions 

The representation of solution in the previous section allows us to obtain 
the coefficients of the series (|13 |) .([14 |) up to any given order. Therefore, if 
k > 1/(2 J +1), we may calculate asymptotic representation of the solution 
with any desired accuracy for r — > 00. Nevertheless, this does not mean that 
series (II Hj) and (fUjl are convergent. Therefore, it would be useful to check 
the results of the previous section using the self-similar solutions, when we 
deal with ordinary differential equations. 

We choose the similarity variable as £ = r/t. In case of the EOS p = 
ke — Ci/(1 + k) the system ©-(0) leads to the ordinary differential equations 

Jv 2 pjjl -v 2 ) +v(l -vQ 

di \{v - if - «(1 - o 2 \ 

^" ( + } [( v -0 2 -«(i-<) 2 ] ' ( } 

where v = u r /u l = (Au 2 — l)/(4w 2 + 1), a = Inp, e = r v p + c\j{\ + k), 
P = u/[J(1 + k)]. 

A qualitative investigation of self-similar relativistic fluid flows with spher- 
ical and cylindrical symmetry can be found in Here we confine ourselves 
to asymptotical solutions of the Eqs. (J2UJ) . (|2Tj) near the boundary of ex- 
panding fluid, where u » 1, v — > 1. If the fluid expands to vacuum, its 
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boundary moves with the light speed (correspondingly, £ = 1). Using the 
results of qualitative invesigation of (3] in the phase plane {f, £}, it is easy 
to classify possible asymptotic solutions of Eqs. (j2*UJl . (j2"Tj) near this boundary 
for £ — > 1 and to obtain their asymptotics. 
In case of n > 1/(2 J + 1) we have 

u 2 ^A 1 (l-C) p \ Pl = -^-, (22) 

e*A 2 r»(l-0 P2 , P2 = v+ Ji \ +K) , (23) 

A\,Ai being arbitrary constants. 

This corresponds to the results of the previous section, where one should 
assign 

Poo(a) ~ a P2 , Voo(a) ~ a u+J . 

If k < 1/(2 J + 1), this type of solution is absent. 

On the other hand, there is also a solution with asymptotics 

u 2 = A ± (l-0'\ A± = (1-k)-\B±Vd), (24) 

where 

B = 1 + k + Jk{1 + 2/i), D = B 2 + (1 - k)[(1 + 2J)k - 1]. 

For k > 1/ (2J + 1) it is easy to see that A_ < 0, so only A + > survives. 
For k < 1/(2 J + 1) both A_ and A + are positive. Depending upon v and 
k, Eq . (|24j) may represent either a special solution or an infinite family of 
solutions. Some values of the parameters may lead to e — > oo for £ — ► 1 
and they must be matched to another solutions through the shock wave; the 
other parameters correspond to solutions that can be matched to vacuum 
(e — > 0). This asymptotics is not well described by a finite number of terms 
in the representation (fT3j) . (JHJ) of the previous section. 

6 Discussion 

We studied non-stationary ultra-relativistic hydrodynamical flows having 
spherical or cylindrical symmetry in case of the equation of state p = p(e). 
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The results may be applicable also in case of a more general EOS if the fluid 
motion is isentropic. The solutions derived in Sect. 4 (as well as the first 
order approximations of Sect. 3) contain two arbitrary functions, therefore 
they represent some kind of a general solution of the equations of the fluid 
motion. However, consideration of the section 5 shows that they do not de- 
scribe solutions of all possible types equally well in the ultra-relativistic limit. 
Formally we may use some of the above results in case of arbitrary k G (0, 1), 
p(e) — ke — Ci, in any ultra-relativistic case, but in certain limited domain 
of spatial variables, where the approximations work. 

The approximations work in the best way in case of the "stiff" equation 
of state: k > 1/5 in case of a spherical expansion ( this involves the im- 
portant 1/3), and k > 1/3 in case of a cylindrical one. In these 
cases the representations of Sect. 3 are most effective for t — r << t — > oo. 
Corresponding solutions exhibit the wave-like behavior irrespectively of en- 
ergy profile of the expanding shell. The energy of the spherical/cylindrical 
shell, which expands almost with the light speed, is approximately conserved. 
This does not mean, however, that this shell will contain all the energy of 
the expanding fluid. It would be interesting to study, what is the part of 
the total energy carried by this expanding wave-like shell, and to relate the 
shell profile with the initial conditions; this probably cannot be done without 
numerical investigation. On the other hand, in case of small values of k we 
may expect that the energy of the fluid is spread over all the values of r ~ t. 
The values of the speed of sound Cg = 1/ (2J+ 1) may be considered as some 
bifurcation points, which separate different regimes of ultra-relativistic flow. 

The method of Section 3 allows to derive formally the solutions up to 
any order of approximation. However, we do not make here any statement 
about the convergence of the series, taking in mind that calculation of high 
orders involves higher order derivatives. Nevertheless, our results show that, 
in spite of some degeneracy in the main order of hydro dynamical equations in 
ultra-relativistic limit, we do not meet any problems like small denominators 
in the higher orders and in this sense we may be sure about the asymptotic 
properties of the solutions. 
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